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A classical nova occurs when material accreting onto the surface of a white dwarf in a close 
binary system ignites in a thermonuclear runawaj^i^l^. Complex structures observed in the 
ejecta at late stageiP® could result from interactions with the companion during the com¬ 
mon envelope phas^^. Alternatively, the explosion could be intrinsically bipolar, resulting 
from a localized ignition on the surface of the white dwarf or as a consequence of rota¬ 
tional distortionPi^. Studying the structure of novae during the earliest phases is challenging 
because of the high spatial resolution needed to measure their small sized^^l Here we report 
near-infrared interferometric measurements of the angular size of Nova Delphini 2013, start¬ 
ing from one day after the explosion and continuing with extensive time coverage during the 
first 43 days. Changes in the apparent expansion rate can be explained by an explosion model 
consisting of an optically thick core surrounded by a diffuse envelope. The optical depth of 
the ejected material changes as it expands. We detect an ellipticity in the light distribution, 
suggesting a prolate or bipolar structure that develops as early as the second day. Com¬ 
bining the angular expansion rate with radial velocity measurements, we derive a geometric 
distance to the nova of 4.54 ± 0.59 kpc from the Sun. 

Nova Delphini 2013 (also known as V339 Del) was discoverecP^I by Koichi Itagaki on 2013 
Aug 14 Universal Time (UT) at 14:01. We began an intensive observing campaign to measure the 
size of the nova with the CHARA Arraj^^, an optical/infrared interferometer located on Mount 
Wilson, California. Our observations began within 15 hours of the discovery and within 24 hours 
of the detonation itself. We measured the size of the expanding ejecta as the nova rose to peak 
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brightness and continued monitoring it for a total of 27 nights between UT 2013 Aug 15 and Sep 
26. 


We measured the angular diameter of Nova Del 2013 by fitting a uniformly bright circle 
to the visibility amplitudes of the interference fringes recorded during each night (Extended Data 
Table 1). We plot the expansion curve in Figure 1. The dotted line shows a linear fit to the 
angular diameters during the first 27 days after the explosion. The inset panel shows an apparent 
deceleration during the first four nights compared with lines of constant velocity. During the last 
week, the measurements show a large jump in the effective size of the nova compared with an 
extrapolation of the linear fit. For comparison, we show the visible and infrared light curves in 
Figure 2. 

The apparent changes in slope of the expansion curve can be explained by a two-component 
model consisting of an optically thick pseudo-photosphere surrounded by an optically thin halo. 
We approximated the intensity distribution projected on the sky using a uniformly bright, circular 
core surrounded by a circular ring. We fit the two-component model to the interferometric visibil¬ 
ities and minimized the total simultaneously across all nights. We fixed the time of detonation 
at modified Julian day to = MID 56518.277, computed by extrapolating the first two pre-discovery 
photometric measurement^ back to the quiescent fluxpEl From our interferometric data, we 
measured an expansion rate of 0.156 ± 0.003 milli-arcseconds (mas) per day for the core diameter 
and a size ratio of 1.73 ± 0.02 between the outer ring and core. We allowed the flux ratio between 
the ring and core to vary on a nightly basis. 

Figure 3 shows how the percentage of light from the ring changes over time. During the 
first two nights, the core and ring have similar surface brightnesses. Therefore, at the earliest 
times, the nova can be approximated by a single uniform disk component where the optically thick 
pseudo-photosphere extends to the outermost, fastest-moving layers of the ejecta. After the peak 
in the visible light curve, the amount of flux in the ring drops, which indicates that the outer layers 
became optically thin and the pseudo-photosphere moved toward the inner, slower-moving, and 
denser layers. Therefore, the apparent deceleration in the expansion curve during the first few 
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nights is likely caused by a diminishing contribution from the outer layers. Over the next few 
weeks, the absolute flux of both components decreased, however, the flux from the core dropped 
at a steeper rate compared with the ring. Therefore, the percentage of light in the ring increased 
relative to the core. These changes are consistent with the spectral evolutiorPI of Nova Del 2013. 
Optical spectra initially showed P Cygni absorption features during the optically thick fireball stage 
which disappeared after about five days. Afterwards, strong emission lines began to dominate the 
spectrum, which are thought to form in the outer, optically thin layerP. The rising strength of the 
emission lines relative to the continuum was also seen in spectra at infrared wavelengths (Extended 
Data Fig. 1). In the two-component model, we assumed a constant ratio between the size of the 
core and the halo. In reality, changes in the optical depth of the ejecta are likely to be more complex 
and could be investigated further by fitting a physical modeP^El simultaneously to interferometric 
and spectroscopic data. Understanding how the optical depth changes within the expanding ejecta 
impacts the physical interpretation of novae light curved. 

During the last five nights, the two-component model shows an enhancement in the flux 
of the outer layers, with the ring contributing ~ 68% of the total light. This coincides with an 
increase in the infrared flux from the novaP^ (see Fig. 2). This suggests the formation of dust in 
the outer layers, perhaps in denser clumps within the ejectaP, that contributes thermal emission. 
Alternatively, as the ejecta expands, a recombination front might propagate through the materiaP^. 
Near-infrared free-free emission from ionized gas the inner regions would decline and shift the 
effective size boundary to the outer layers. This interpretation is consistent with the larger size of 
the effective photosphere at later dates, however, it does not explain the brightening in the near- 
infrared light curve. The effective size of our two-component model is shown as the solid line in 
Figure 1; the model reproduces the changes in slope of the expansion curve. 

We reconstructed model-independent images of Nova Del 2013 using the MArkov Chain 
IMageil^ (MACIM) for data obtained on three nights during the first week with sufficient baseline 
coverage. As shown in Figure 4, the images show a striking similarity to the two-component 
model, with an optically thick core surrounded by a halo of fainter emission. 
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Combining the angular expansion rate with the radial velocity of the ejected material pro¬ 
vides a way to measure the distanced to Nova Del 2013. The velocities reported in the literature!^ 
range from —600 to —2500 km s“^, consistent with the expectation that lines form at different 
layers in the expanding atmosphere with a radial outflow velocity proportional to the distance from 
the white dwarf explosion site. We selected the Si II A A 6347, 6371 A absorption lines observed 
during the first week to represent the deeper layers, because they likely form in the dense, cooler 
gas immediately above the continuum forming layer. Based on an analysis of spectra downloaded 
from the archive of the Astronomical Ring for Access to Spectroscop)^^, we estimated the outflow 
speed near the continuum forming layer of Kjection = 613 ± 79 km s“^ Combining this veloc¬ 
ity with the angular expansion rate, we derived a distance of 4.54 ± 0.59 kpc to Nova Del 2013. 
Empirical evidence suggests that novae have prolate spheroid^ or bipolai^SI shapes, therefore the 
true distance will depend on its geometry and inclination in the plane of the sky!^. Assuming a 
probable inclination of 45° for Nova Del 2013 (S. Shore, personal communication, 16 April 2014), 
we estimate that these effects contribute to a systematic difference in the distance of less than 
5%. Using the light curve decay times and reddening reporter?^ for Nova Del 2013 (assuming an 
Ry = 3.1 extinction law), the maximum magnitude - rate of decline (MMRD) relatiorP^ gives a 
distance ranging from 3.3 to 4.1 kpc, depending on the adoption of the linear or non-linear MMRD 
formulations. The expansion parallax also confirms the distance of 4.2 kpc based on a spectral 
comparisorP^to nova OS Andromedae 1986. 

To investigate the development of asymmetric features in Nova Del 2013, we fit a uniformly 
bright ellipse to the visibilities on nights with coverage over a range of position angles on the sky 
(from two to nine days after the outburst; Extended Data Table 2). Eor each night, the ellipse 
produces a significantly lower compared with the circular disk (P < 0.01 for most nights; 
see Methods). The nova appears to be ~ 13% larger along the major axis compared with the 
minor axis, suggesting that the ejecta could be bipolar. Eor the remaining nights with limited sky 
coverage, we fixed the axis ratio and position angle to their average values (6*major/^miiior = 1.13 ± 
0.07, PA = 128° ± 30°) and solved for the size of the ellipse. On average, we find that the circular 
disk diameters underpredict the mean angular diameter of the ellipse by ~ 2%. Therefore, the 
assumption of circular symmetry in the two-component model should not strongly bias our results. 
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In addition to the visibility amplitudes that provide information on the size and shape of 
the source, interferometric closure phases obtained when combining the light from three or more 
telescopes indicate deviations from point symmetry (whether the brightness is the same when 
reflected through a point at the center of the source distribution). The closure phases from the 
first two sets of data obtained with all six telescopes (t = 3.0 and 5.0 days) are consistent with 
zero within the uncertainties, suggesting that the source was point-symmetric within the resolution 
limits. By the third set of observations with all six telescopes (t = 6.9 days; see Fig. 4), the closure 
phases rise steadily to ~ 60° (Extended Data Fig. 2), suggesting the detection of a point asymmetry 
in the brightness distribution. If the outflow is bipolar, as suggested by our elliptical fits, then this 
could indicate a difference in the brightness between the two lobes (perhaps due to viewing angle) 
or, alternatively, the development of clumpy structures within the expanding material. 

The ellipticity and closure phase asymmetries indicate that non-spherical structures devel¬ 
oped as early as a few days after the outburst. Along with results from previous interferometric 
studies of novaeP®^, this suggests that novae explosions might be inherently bipolar or that the 
elliptical shape develops early during the common envelope phase. 
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Full Methods and any associated references are available in the online version of the paper at 
www.nature.com/nature. 
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Figure 1 - Expansion curve of Nova Del 2013. The angular diameters were measured by fitting 
a circular disk to the interferometric data. The dotted line shows a linear fit for days 0-27. The 
inset panel zooms in during the first week and shows dotted lines with velocities of 613 and 2500 
km s“^ at a distance of 4.54 kpe. The apparent deeeleration during the first week, along with the 
jump in size during the last week, ean be explained by a two-component model consisting of a 
circular core surrounded by a ring where the flux ratio changes over time. The grey region shows 
the expansion rate of the core (lower edge) and the ring (upper edge). The solid line shows the 
effective size obtained by fitting the two-component model visibilities as a single circular disk. 
Error bars represent 1 cr uncertainties computed from a boostrap analysis. 
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Figure 2 - Infrared light curve of Nova Del 2013. The magnitudes and error bars were computed 
from the mean and standard deviation of counts recorded on the detector during each interferomet¬ 
ric observation collected at the CHARA Array. The data were obtained in the i7-band using two 
different beam combiners (blue squares and green diamonds) and in the /f'-band (red triangles). 
For comparison, we plot daily averages of photometric measurements in the V -band downloaded 
from the American Association of Variable Star Observers (AAVSO; black circles). The rise to 
peak brightness was slower in the infrared compared to the visible. 
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Figure 3 - Changes in the flux ratio of the two-component model. The solid line shows a fit 
where the ratio changes smoothly over time. For the first three nights, the ring represents the 
outer boundary of the optically thick pseudo-photosphere. After the peak in the visible light curve, 
the flux in the ring drops. On nights 4-27, the percentage of light from the ring increases at an 
approximately linear rate. On the last four nights, the ring contributes an average of 68% of the 
total light. Error bars represent 1 a uncertainties derived from a least squares fit to the visibility 
data. 
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Figure 4 - Model and reconstructed images of Nova Del 2013. The two-component model consists 
of a circular core surrounded by a fainter ring (top row). The images reconstructed using MACIM 
show a similar spatial structure (bottom row). The data used in the reconstructions were obtained 
with the Michigan Infrared Combiner (MIRC) at the CHARA Array on UT 2013 Aug 17, 19, and 
21 (t = 3.0, 5.0, and 6.9 days after the outburst) and have good baseline coverage on the sky. We 
scaled the flux by the infrared i7-band magnitude measured on each night to show how the surface 
brightness changes. 
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Methods 


Interferometric Observations with the CHARA Array 

The CHARA Arrajff^ is an optical/infrared interferometer located on Mount Wilson, Cal¬ 
ifornia. The array has six 1 m telescopes arranged in a Y configuration with baselines ranging 
from 34 to 331 m. The longest baseline provides a spatial resolution of 0.5 mas at a wavelength of 
1.6 /im. We observed Nova Del 2013 in the infrared using the CLASSIC two-beam combiner^\ 
the CLIMB three-beam combiner^^, and the MIRC six-beam combiner^^. The data were reduced 
using the standard reduction pipeline for each instrument. Data products include the squared, 
normalized visibility amplitude of the interference fringes for the two-baseline combiners and 
the squared,normalized visibility amplitudes and closure phases for the instruments that combine 
the light from three or more telescopes (see Extended Data Figures 2-4). Before and after each 
observation, we observed unresolved stars or stars with known angular diameters^^’^^’^® to cali¬ 
brate the interferometric measurements. Extended Data Table 1 lists the UT date, modified Julian 
Day (MID), time since the outburst (to = MID 56518.277), beam combiner, telescope configu¬ 
ration, filter, and calibrators used during each observation. The calibrated data will be available 
though the Optical Interferometry Database developed by the Jean-Marie Mariotti Center (JMMC; 
http://www.jmmc.fr/oidb.htm). On five nights we also collected data using two visible light beam 
combiners; these data will be discussed in a subsequent publication (O.C., D.M., RT., V.M., I.T-B, 
D.P.K.B, E.Lagadec, M.I., N.N., A.M, P.S, R.T.Z., G.H.S, T.t.B., D.R.G., H.A.M., T.B., C.D.F., 
N.V., N.S., J.S., and L.S., manuscript in preparation). 

In addition to the interferometric measurements, we also computed infrared magnitudes of 
Nova Del 2013 using the number of counts recorded on the detector. We compared the counts from 
Nova Del 2013 with the calibrator stars observed immediately before and after. We calibrated the 
photometry based on the 2MASS magnitudes^® of the calibrators. We present the photometry in 
Extended Data Table 3. Figure 2 shows the infrared light curve of Nova Del 2013. For comparison, 
we also plot daily averages of C-band photometric measurements downloaded from AAVSO®^. 
The rise to peak brightness was slower in the infrared compared to the visible. Based on a second- 
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order polynomial fit to the H and it^-band measurements near the maximum (< 4.3 mag), the 
near-infrared light curve reached a maximum of 3.61 mag at MJD 56524.6 (t = 6.3 days). The 
maximum brightness in the visibleP^ occurred at an earlier time of MJD 56520.4 (t = 2.1 days) at 
V ~ 4.46 mag. 


Correcting for the Effects of Line Emission 


The emergence of strong emission lines in the infrared spectrum of Nova Del 2013 (such as 
Br 7 in the /f-band) will broaden the central envelope of the interference fringe and create sidelobes 
farther out^®. We investigated the effects of the emission lines on the calibration of the interfero¬ 
metric data obtained with the CLASSIC and CLIMB beam combiners. The reduction code outputs 
the mean power spectrum for each target. We fit a rectangular top-hat function to each power spec¬ 
trum and calculated the mean center and width for the observations on each night. By comparing 
the results from the nova with the calibrators, we found that the central wavelength of the filter 
remained constant over time. We also found that the effective bandwidth of the filter remained the 
same during the first two weeks and then decreased by ~ 14% ± 3% as the emission lines ap¬ 
peared in the nova spectra (see Extended Data Table 4). Essentially, the emission lines put a larger 
fraction of the flux at the center of the wavelength region, thereby effectively decreasing the full- 
width at half-maximum (EWHM) of the filter. The CLASSIC/CLIMB reduction software assumes 
a constant top-hat function for the width of the filter. Therefore, when the effective bandwidth 
decreases, the visibility is overestimated and the angular size of the nova is underestimated. 

To correct for this effect, we multiplied the CLASSIC/CLIMB visibilities by the ratio of 
the mean effective width of the nova observations relative to the calibrators measured on each 
night. We added the uncertainty in the bandwidth ratio in quadrature with the uncertainties in the 
visibilities reported by the reduction code. We applied this correction starting on UT 2013 Aug 
28, when the effects of the emission lines became measurable. During the first two weeks, when 
the contribution from the emission lines had a negligible effect on the calibration, we assumed the 
default bandwidth used in the reduction code, but included the average uncertainty in the nova-to- 
calibrator bandwidth ratio of ±0.164 in the error budget so that all of the CLASSIC and CLIMB 
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visibilities would have similar weights during the global two-component fit. For the uniform disk 
fits to the data from each night, the bandwidth correction led to the measurement of the angular 
size of the nova being larger at the ~ 5% level after the emission lines developed. 


Uniform Disk and Uniform Ellipse Fits 


For each individual night of CHARA data, we measured the angular diameter of Nova Del 
2013 by fitting a uniformly bright, circular disk to the visibilities and performing a Levenberg- 
Marquardt least-squares minimization using the IDL mpfit package^®. The uniform disk diameters 
(0 ud)> reduced xt (where u is the number of degrees of freedom), and number of visibility mea¬ 
surements N{V‘^) on each night are listed in the last three columns of Extended Data Table 1. 

For five nights with sufficient baseline coverage over a range of position angles, we fit a 
uniformly bright ellipse to the visibility data. The fitted parameters include the size of the major 
and minor axes and position angle of the major axis (6*major, (*minor> PA); the values are given in 
Extended Data Table 2. For four out of the five nights, performing an F-test indicates that the 
uniform ellipse produces a significant improvement in the x^ compared with a uniform circle, at a 
significance level of 0.01. On UT 2013 Aug 16, the improvement is only at the 0.10 significance 
level. 


For both the uniform disk and uniform ellipse fits, we determined uncertainties in the pa¬ 
rameters using a bootstrap approach^®. During each iteration, we created a synthetic sample of 
visibility measurements from our original set of measurements by randomly selecting the same 
total number of visibilities, with replacement. Therefore, in the synthetic sample, some measure¬ 
ments are repeated and others are left out. We applied normally distributed uncertainties to the 
synthetic measurements (using the measured uncertainties from the original data). We then fit a 
uniform circle or uniform ellipse to the synthetic data set to determine the best fit parameters for 
that iteration. We performed 1000 iterations and computed uncertainties based on the standard de¬ 
viation of the parameter distributions. These uncertainties are reported in Extended Data Tables 1 
and 2. 
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We expected the size of the nova to vary over the course of the night. The maximum length 
of observing time was ~ 6 hours on UT 2013 Aug 16 and Sep 6. Given an expansion rate of 
0.156 mas/day, we would expect the size of the nova to change by ~ 0.04 mas over the course 
of the observations. For the three nights with the smallest uncertainties on the angular diameter 
measurements (± 0.005 mas on UT 2013 Aug 17,19, and 21), the length of observing time ranged 
from 30 to 45 minutes, so the change in size is expected to be < 0.005 mas. For all other nights, 
the length of observing time varied from ~ 30 minutes to 3 hours, corresponding to a change in 
size of 0.003 to 0.02 mas. For these nights, the quoted uncertainties on the angular diameters in 
Extended Data Table 1 are larger than ± 0.02 mas, which should account for any changes in size 
over the course of the night. For the two-component model described in the next section, we fit a 
dynamic model directly to the visibilities versus time, so it accounts for changes in size over the 
course of the night. 


Two-Component Model 


The two-component model consists of a uniformly bright, circular disk surrounded by a uni¬ 
formly bright ring. The angular diameter of the central uniform disk (^ud) and the outer diameter 
of the ring (^ring) are given by: 

(if) 

= ( 1 ) 

^rmg(f) = Coring ^UD(f) (2) 

where ^ is the angular expansion rate of the central core and is the ratio of the outer ring 
diameter relative to the uniform disk diameter. 

We determined the two global parameters and Cring) by minimizing the cumulative 
between the measured and model visibilities across all nights. To do this, we searched through 
a grid of values for the expansion rate and size ratio and calculated the size of each component 
during the times of observation. Then, for each night of observation, we performed a Levenberg- 
Marquardt least-squares minimization using the IDL mpfit package^® to solve for the best-fitting 
flux ratio between the ring and the core component. 
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There are degeneracies between the overall size of the nova and the flux ratio of the two 
components; a smaller expansion rate can be compensated by placing more flux into the outer 
component. Therefore, we followed a two-step procedure in fitting the expanding two-component 
model simultaneously to all nights. Initially, we allowed the flux in the ring relative to the core 
to vary each night, but with a limit set so that the surface brightness of the ring would not exceed 
that of the core. The initial restriction prevented a large increase in the fitted brightness of the ring 
at later nights. Using the grid-search described above, we determined an expansion rate of 0.156 
± 0.003 mas per day for the core diameter and a size ratio of 1.73 ± 0.02 between the outer ring 
diameter and the core diameter. We then fixed the expansion rate and size ratio and solved for the 
ring-to-core flux ratio for each night after removing the restriction on the surface brightness on all 
but the first three nights (we would not expect a clearing of the central region at such early times). 
This allows for a moderate brightening of the ring at later times. 

In Figure 1, the solid line shows the effective angular size of the two-component model. 
This curve was generated by assuming that the percentage of light from the ring varies smoothly 
over time (solid line in Fig. 3) and computing the model visibilities for each night based on the 
expansion rate and size ratio. We then fit a single uniform, circular disk to the model visibilities 
from each night to compute the effective size. Extended Data Figure 5 shows how the size and flux 
ratios of the two-components change during the nights of observation. 

Emission Line Flux from Infrared Spectroscopy 

To investigate how the emission line flux changed in the infrared, we used low dispersion 
spectra obtained with the TripleSPEC spectrograph"^^’^^ on the 200-inch Hale Telescope at Palomar 
Observatory on UT 2013 Aug 20 and 23. We reduced the spectra using the Spextool reduction 
program modified for TripleSpec^^, which includes routines to correct for detector artifacts^^. We 
calibrated and removed the telluric lines in the nova spectrum using the xtellcor program^^, which 
compares a high signal-to-noise spectrum of Vega to observations of a nearby AO star, SAO 88500, 
taken immediately after each observation of the nova. We also used infrared spectra^® obtained at 
Mount Abu, India^^ on UT 2013 August 28 and 29, and September 8 and 20. 
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We collected these spectra into two regions covering the wavelength range of the CHARA 
CLASSIC/CLIMB i/-band and A''-band filters and normalized the flux to unity at 1.65 and 2.2 
/rm, respectively. The i7-band is dominated by the upper level transitions of the hydrogen Brackett 
series plus a number of C I lines, and the A''-band is dominated by Br 7 , He I A 2.0585 /um, and 
some C I lines^^. 

To estimate the flux from the continuum, we made spline fits to the lower envelope of each 
spectrum. We then summed both the continuum flux and the emission flux (total minus continuum) 
across the wavelength bands to form a ratio of the integrated emission line flux relative to the 
integrated continuum flux. As shown in Extended Data Figure 1, the emission strength grew over 
the time-frame of the CHARA observations of Nova Del 2013 and began to decrease in the final 
spectrum. These changes are similar to the changes we measured in the effective bandwidth of the 
interferometric observations. 

The strengthening in the emission lines helps explain why the expansion curve of Nova Del 
2013 became steeper after the first week. The emission lines form over a greater radial extent 
in the expanding envelope than does the continuum from the pseudo-photosphere. Thus, if the 
outflow follows a Hubble-type relation with the velocity proportional to radial distance^®, then the 
emission lines sample parts of the envelope that were expanding faster than the photosphere. We 
interpolated linearly between the spectroscopic measurements to obtain the ratio of the emission 
line to continuum flux at the times of the CHARA observations and used this as the initial estimate 
of the flux ratio in the halo ring relative to the uniform disk core in the two-component fit. 

The drop in the emission to continuum flux ratio after the fourth week can be attributed to 
the onset of dust formation. This can be seen in the A'-band spectra, where the continuum slopes 
down in the first spectra and becomes flat in the last spectrum. This is consistent with flux from 
a dusty component that contributes more to the continuum and increases with wavelength. Our 
two-component model based on the interferometric observations suggests that if dust formation 
occurred, that it happened over a larger spatial scale than that of the expanding photosphere. 

We also made measurements of the flux weighted central wavelength of each infrared spec- 
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trum. These were essentially constant to within a percent or so during these observations, consis¬ 
tent with our results from the analysis of the fringe power spectra. 

Measuring the Radial Velocity of the Optically Thick Core 

In order to compare the angular expansion rate with the radial velocity of the nova outflow, 
we needed to decide which features in the optical spectra are most representative of the kinemat¬ 
ics of the near-infrared continuum. The velocities reported in the literaturepi range from —600 to 
—2500 km s“^, consistent with the expectation that lines form at different layers in the expanding 
atmosphere with a radial outflow velocity proportional to distance from the white dwarf explosion 
site. The optical depth unity boundary for the near-infrared continuum must form at higher densi¬ 
ties deep inside the outflow (compared to the P Cygni lines and emission lines that form further out 
in lower density gas). The absorption lines observed during the first week are probably the most 
representative of the deeper layers, because they form in the dense, cooler gas immediately above 
the continuum forming layer. We selected the Si II AA 6347, 6371 A absorption lines for measure¬ 
ment because they were observed throughout the first week and are relatively free of blending and 
interference from Earth’s atmospheric lines and nova emission lines. 

We downloaded six spectra with a high resolving power {R - 10,000) obtained over the first 
week by Olivier Garde from the archive of the Astronomical Ring for Access to SpectroscopjP 
and transformed them to a unit continuum on a heliocentric wavelength grid. We measured radial 
velocities for the Si II lines by cross-correlation with a model spectrum (Teff = 8000 K, log g — 2.0, 
solar metallicity) based upon ATLAS atmospheres from Robert Kurucz. The mean velocity and 
standard deviation from this sample (MJD 56519.8 - 56524.8) are —598 ± 51 km s“^ The local 
standard of rest radial velocity at the probable distance of the nova is -1-15 km s“^, so the implied 
outflow velocity is —613 ± 51 km s“^. However, we do not know the peculiar velocity of the 
nova relative to its local standard of rest, so we applied a representative velocity dispersion of 60 
km s“^ for white dwarfs to revise the uncertainty in our final estimate of the outflow speed near the 
continuum forming layer, ^ejection = 613±79 km s~^. Note that this is strictly an upper limit for the 
velocity of the continuum forming layers because the absorption lines form above the continuum 
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in a somewhat faster outflow. Additionally, the estimate reflects the geometry of the outflow which 
is dependent on the axial inclination for a bipolar ejection. For example, if the absorption lines 
formed in a spherically expanding photosphere, then we might expect to observe line boundaries 
ranging from zero (flux from the limb) to — Kutflow (flux from the center), with a mean velocity 
between these limits. However, the observed Si II lines are rather narrow (FWHM ~ 280 km s“^), 
contrary to this prediction. On the other hand, for a bipolar outflow the absorption velocities may 
be much more restricted to the radial ejection speed along the line of sight. Consequently, until 
future observations reveal details of the geometry of the nova remnant, we will simply assume that 
the measured radial velocity is the same as the radial outflow speed in our direction and that this is 
also the same as the transverse expansion speed as observed in the plane of the sky. 
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Extended Data Table 1 - Journal of observations and angular diameters measurements of Nova 
Del 2013. The angular diameters were computed by fitting a circular uniform disk to the visibilities 
from each night. 


UT Date 

MJD 

t (days) 

Comb* 

Configuration 

Filtert 

Ca|t 

OuD ( 

mas) 

xi 

N(V2) 

2013 Aug 15 

56519.231 

1.0 

CL 

S1E2W1 

H, K' 

1,2 

0.414 


0.076 

1.19 

18 

2013 Aug 16 

56520.303 

2.0 

CL 

SlEl 

H 

1,2 

0.537 


0.028 

0.23 

9 

2013 Aug 17 

56521.260 

3.0 

M 

W1S2S1E1E2W2 

H 

3 

0.652 


0.005 

1.40 

180 

2013 Aug 18 

56522.284 

4.0 

c 

SlEl 

H 

1 

0.777 


0.038 

1.77 

4 

2013 Aug 19 

56523.248 

5.0 

M 

W1S2S1E1E2W2 

H 

4,3 

0.863 


0.004 

4.10 

180 

2013 Aug 20 

56524.325 

6.0 

CL 

E2W2W1 

H 

1 

1.096 


0.056 

0.82 

6 

2013 Aug 21 

56525.220 

6.9 

M 

W1S2S1E1E2W2 

H 

4,3 

1.250 


0.006 

11.5 

240 

2013 Aug 22 

56526.216 

7.9 

C 

ElWl 

K' 

4 

1.262 


0.046 

0.17 

2 

2013 Aug 23 

56527.310 

9.0 

CL 

SIEIWI, S2E2W2 

H, K' 

1,4,5 

1.643 


0.077 

5.10 

15 

2013 Aug 24§ 

56528.242 

10.0 

M 

W1S2S1E1E2W2 

H 

3 

1.885 


0.016 

4.86 

32 

2013 Aug 26§ 

56530.208 

11.9 

M 

W1S2S1E1E2W2 

H 

3 

2.314 


0.037 

1.09 

61 

2013 Aug 28 

56532.191 

13.9 

CL 

W1W2 

H, K' 

4 

2.62 


0.27 

2.49 

3 

2013 Aug 29 

56533.235 

15.0 

CL 

W1W2 

K' 

4 

3.133 


0.076 

0.29 

2 

2013 Aug 30 

56534.273 

16.0 

CL 

E1E2, S1S2 

K' 

1,4 

3.89 


0.72 

0.19 

3 

2013 Sep 02 

56537.146 

18.9 

C 

E1E2 

K' 

4 

3.51 


0.31 

1.10 

2 

2013 Sep 04 

56539.178 

20.9 

C 

E1E2 

K' 

1,4 

4.22 


0.22 

3.32 

3 

2013 Sep 05 

56540.164 

21.9 

C 

E1E2 

K' 

1,4 

4.57 


0.18 

0.04 

3 

2013 Sep 06 

56541.280 

23.0 

C 

E1E2 

K' 

1,4 

4.62 


0.13 

0.40 

16 

2013 Sep 07 

56542.158 

23.9 

C 

E1E2 

K' 

1 

4.94 


0.11 

0.23 

6 

2013 Sep 08 

56543.150 

24.9 

CL 

E1E2 

K' 

1 

5.11 


0.12 

1.68 

5 

2013 Sep 09 

56544.156 

25.9 

CL 

E1E2 

K' 

1 

5.49 


0.11 

1.76 

6 

2013 Sep 10 

56545.139 

26.9 

CL 

E1E2 

K' 

1 

6.09 


0.17 

0.38 

2 

2013 Sep 13 

56548.160 

29.9 

CL 

S1S2 

K' 

1 

5.19 


0.97 

0.31 

3 

2013 Sep 16 

56551.137 

32.9 

C 

S1S2 

K' 

1 

9.12 


0.74 

0.85 

3 

2013 Sep 19 

56554.138 

35.9 

C 

S1S2 

K' 

1 

9.33 


0.30 

0.12 

5 

2013 Sep 23 

56558.121 

39.8 

C 

S1S2 

K' 

1 

10.96 


0.18 

0.52 

5 

2013 Sep 26 

56561.159 

42.9 

C 

S1S2 

K' 

1 

11.83 


0.98 

0.01 

2 


Footnotes to Extended Data Table 1 - 

* Beam Combiner Codes: C = CLASSIC, CL = CLIMB, M = MIRC 

t CLASSIC/CLIMB H: X = 1.6731 ^m, AA = 0.2854 ^m; CLASSIC/CLIMB A': A = 2.1329 
fim, AA = 0.3489 fim; MIRC H: X- L5-L7 fim with AA ~ 0.034 yum for each of the 8 spectral 
channels. 

* Calibrator Codes: 1 - HD 190993 (6 = 0.161 ± 0.011 mas [33]), 2 - HD 196740 (6 = 0.180 ± 
0.012 mas [33]), 3 - HD 206860 (6 = 0.530 ± 0.015 mas [34]), 4 - HD 190406 (6 = 0.584 ± 
0.010 mas [35]), 5 - HD 195810 (9 = 0.301 ± 0.021 mas [33]) 
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§ Only E2W1W2 and S1S2 were recorded on UT 2013 Aug 24. Only W1W2, E1E2, S1S2, and 
E2W2 were recorded on UT 2013 Aug 26. 


Extended Data Table 2 - Uniform ellipse models of Nova Del 2013. 


UT Date 

MJD 

t (days) 

Comb 

^major (mas) 

^minor (mas) 

PA (°) 

xl 

2013 Aug 16 

56520.303 

2.0 

CL 

0.631 ± 0.282 

0.503 ± 0.243 

81.8 ± 31.9 

0.14 

2013 Aug 17 

56521.260 

3.0 

M 

0.676 ± 0.011 

0.633 ± 0.007 

114.8 ± 10.6 

1.30 

2013 Aug 19 

56523.248 

5.0 

M 

0.936 ± 0.005 

0.826 ± 0.004 

148.7 ± 2.1 

1.53 

2013 Aug 21 

56525.220 

6.9 

M 

1.291 ± 0.007 

1.198 ± 0.005 

157.4 ± 3.5 

6.95 

2013 Aug 23 

56527.310 

9.0 

CL 

1.716 ± 0.094 

1.535 ± 0.044 

138.5 ± 16.1 

2.49 
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Extended Data Table 3 - Infrared Magnitudes of Nova Del 2013. 


MJD 

Filt 

Magnitude 

MJD 

Filt 

Magnitude 

MJD 

Filt 

Magnitude 

56519.179 

H 

5.12 


0.16 

56532.170 

K' 

4.44 


0.23 

56542.149 

K' 

5.37 


0.16 

56519.193 

H 

5.19 


0.17 

56532.191 

K' 

4.53 


0.16 

56542.158 

K' 

5.37 


0.16 

56519.216 

H 

5.08 


0.15 

56532.214 

H 

4.65 


0.38 

56542.167 

K' 

5.35 


0.15 

56519.231 

H 

5.20 


0.14 

56533.223 

K' 

4.65 


0.21 

56542.176 

K' 

5.32 


0.16 

56519.250 

K' 

4.97 


0.17 

56533.235 

K' 

4.60 


0.20 

56543.125 

K' 

5.42 


0.22 

56519.261 

K' 

5.04 


0.17 

56534.249 

K' 

4.74 


0.16 

56543.135 

K' 

5.44 


0.20 

56519.273 

K' 

4.93 


0.16 

56534.265 

K' 

4.78 


0.17 

56543.150 

K' 

5.42 


0.19 

56520.187 

H 

4.20 


0.20 

56534.273 

K' 

4.71 


0.16 

56543.162 

K' 

5.36 


0.19 

56520.218 

H 

4.24 


0.19 

56534.404 

K' 

4.49 


0.27 

56543.175 

K' 

5.30 


0.18 

56520.247 

H 

4.16 


0.17 

56537.130 

K' 

5.08 


0.17 

56544.123 

K' 

5.32 


0.22 

56520.278 

H 

4.21 


0.16 

56537.146 

K' 

4.79 


0.18 

56544.132 

K' 

5.44 


0.19 

56520.303 

H 

4.15 


0.15 

56537.157 

K' 

4.46 


0.31 

56544.145 

K' 

5.52 


0.19 

56520.349 

H 

4.08 


0.16 

56539.168 

K' 

5.10 


0.19 

56544.156 

K' 

5.52 


0.17 

56520.375 

H 

4.11 


0.16 

56539.179 

K' 

5.19 


0.22 

56544.170 

K' 

5.42 


0.15 

56520.412 

H 

4.07 


0.18 

56539.190 

K' 

5.62 


0.23 

56544.184 

K' 

5.37 


0.17 

56520.432 

H 

4.04 


0.19 

56540.148 

K' 

5.25 


0.20 

56545.128 

K' 

5.25 


0.29 

56521.241 

H 

4.09 


0.12 

56540.164 

K' 

5.22 


0.18 

56545.139 

K' 

5.46 


0.30 

56521.263 

H 

4.05 


0.12 

56540.179 

K' 

5.17 


0.17 

56545.153 

K' 

6.56 


0.54 

56522.265 

H 

3.74 


0.15 

56540.194 

K' 

5.14 


0.18 

56548.144 

K' 

5.38 


0.21 

56522.275 

H 

3.82 


0.14 

56540.212 

K' 

5.18 


0.16 

56548.160 

K' 

5.26 


0.21 

56522.284 

H 

3.75 


0.14 

56541.143 

K' 

5.29 


0.15 

56548.171 

K' 

5.09 


0.22 

56522.295 

H 

3.70 


0.14 

56541.160 

K' 

5.31 


0.15 

56551.119 

K' 

5.43 


0.27 

56523.219 

H 

3.63 


0.15 

56541.178 

K' 

5.34 


0.15 

56551.137 

K' 

5.49 


0.26 

56523.248 

H 

3.46 


0.16 

56541.195 

K' 

5.31 


0.14 

56551.154 

K' 

5.28 


0.22 

56524.295 

H 

3.46 


0.16 

56541.213 

K' 

5.22 


0.14 

56554.118 

K' 

5.00 


0.25 

56524.325 

H 

3.52 


0.23 

56541.231 

K' 

5.30 


0.12 

56554.129 

K' 

5.08 


0.23 

56525.189 

H 

3.79 


0.13 

56541.247 

K' 

5.29 


0.13 

56554.139 

K' 

5.08 


0.22 

56525.220 

H 

3.80 


0.14 

56541.265 

K' 

5.34 


0.13 

56554.150 

K' 

5.11 


0.22 

56526.207 

K' 

3.58 


0.15 

56541.281 

K' 

5.30 


0.12 

56554.159 

K' 

5.17 


0.22 

56526.224 

K' 

3.64 


0.13 

56541.299 

K' 

5.22 


0.14 

56558.102 

K' 

4.77 


0.22 

56527.258 

K' 

3.86 


0.13 

56541.315 

K' 

5.10 


0.14 

56558.112 

K' 

4.80 


0.20 

56527.283 

K' 

3.68 


0.13 

56541.333 

K' 

5.41 


0.12 

56558.121 

K' 

4.77 


0.20 

56527.310 

K' 

3.84 


0.16 

56541.350 

K' 

5.19 


0.13 

56558.131 

K' 

4.79 


0.20 

56527.324 

K' 

3.67 


0.16 

56541.367 

K' 

4.86 


0.15 

56558.139 

K' 

4.81 


0.19 

56527.394 

H 

4.29 


0.15 

56541.389 

K' 

5.20 


0.14 

56561.129 

K' 

4.87 


0.35 

56528.241 

H 

3.91 


0.16 

56541.409 

K' 

4.49 


0.17 

56561.159 

K' 

4.54 


0.28 

56530.203 

H 

4.35 


0.12 

56542.126 

K' 

5.38 


0.22 






56530.263 

H 

4.54 


0.11 

56542.141 

K' 

5.29 


0.17 
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Extended Data Table 4 - Effective bandwidth of the CLASSIC/CLIMB observations. We present 
the ratio of the effective bandwidth of the interferometric observations of Nova Del 2013 compared 
with the calibrator stars. 


UT Date 

MJD 

t (days) 

Bandwidth (Nova) 
Bandwidth(Cal) 

2013 Aug 15 

56519.231 

1.0 

0.916 

± 0.295 

2013 Aug 16 

56520.303 

2.0 

1.015 

± 0.256 

2013 Aug 22 

56526.216 

7.9 

0.970 

± 0.193 

2013 Aug 28 

56532.191 

13.9 

1.023 

± 0.134 

2013 Aug 29 

56533.235 

15.0 

0.905 

± 0.107 

2013 Aug 30 

56534.273 

16.0 

0.879 

± 0.151 

2013 Sep 02 

56537.146 

18.9 

0.863 

± 0.104 

2013 Sep 04 

56539.178 

20.9 

0.844 

± 0.041 

2013 Sep 05 

56540.164 

21.9 

0.844 

± 0.121 

2013 Sep 06 

56541.280 

23.0 

0.876 

± 0.208 

2013 Sep 07 

56542.158 

23.9 

0.879 

± 0.117 

2013 Sep 08 

56543.150 

24.9 

0.817 

± 0.078 

2013 Sep 09 

56544.156 

25.9 

0.878 

± 0.086 

2013 Sep 10 

56545.139 

26.9 

0.792 

± 0.166 

2013 Sep 13 

56548.160 

29.9 

0.836 

± 0.098 

2013 Sep 16 

56551.137 

32.9 

0.933 

± 0.194 

2013 Sep 19 

56554.138 

35.9 

0.883 

± 0.162 

2013 Sep 23 

56558.121 

39.8 

0.902 

± 0.111 

2013 Sep 26 

56561.159 

42.9 

1.165 

± 0.492 
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Extended Data Figure 1 - Integrated emission line to continuum flux ratios measured from in¬ 
frared spectroscopy. The squares represent the i/-band ratios while the diamonds represent the 
it''-band ratios. The rise in the emission line flux is consistent with an increasing contribution 
from optically thin emission. The down turn in the curve during the last measurement is likely 
caused by a rising contribution of the continuum from the formation of dust. 
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Extended Data Figure 2 - Closure phases measured with MIRC on UT 2013 Aug 21. The mea¬ 
surements are plotted against the length of the maximum baseline {B) for the group of three tele¬ 
scopes that were used to form the closure phase. The non-zero closure phases indicate an asym¬ 
metry in the light distribution that is not point symmetric. The colors of the symbols are used 
to differentiate the measurements from each grouping of three telescopes in the array. The clo¬ 
sure phases are measured in eight wavelength channels. Error bars represent 1 a measurement 
uncertainties. 
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Extended Data Figure 3 - Example of MIRC observations obtained on UT 2013 Aug 21. The 
light is dispersed over 8 wavelength channels in the iif-band. a, Coverage of interferometric base¬ 
lines projected on the plane of the sky in right ascension (RA) and declination (Dec.) in units of 
spatial frequencies {u = B^/K v = By/X). b. Squared, normalized visibility amplitude measure¬ 
ments, color coded to match the baselines on the left. The solid line shows the best fit uniform 
disk model. The small black dots show the best fit uniform ellipse model. Error bars represent 1 a 
measurement uncertainties. 


30 









Extended Data Figure 4 - Interferometric visibilities of Nova Del 2013 measured with the CHARA 
Array. The red line shows the best-fitting model for a uniformly bright, circular disk. The time 
since the explosion (in days) is indicated in the upper right comer of each panel. The measurements 
were obtained with the CLASSIC, CLIMB, and MIRC beam combiners (see Extended Table 1). 
Error bars represent 1 cr measurement uncertainties. 
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Extended Data Figure 5 - Time evolution of the two-component model of Nova Del 2013. The 
model consists of a circular core surrounded by a halo ring. The expansion rate and size ratio were 
determined by minimizing across all of the nights while allowing the flux ratio to vary night by 
night. In plotting the images, we used the flux ratio measured directly for the first three nights, the 
linear fit plotted in Figure 3 of the main paper for t = 4-27 days, and our measurement that the ring 
contributes an average of 68% of the light on the last five nights (dust emission in the outer layers). 
Each panel is 12 mas on a side. We scaled the model flux by the infrared magnitude measured on 
each night to show how the surface brightness changes. The time since the explosion (in days) is 
indicated in each panel. Intensity refers to the flux per unit area. 



